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ABSTRACT 

We present a photometric search for objects with point-source components that are optically variable 
on timescales of weeks-months in the Hubble Ultra Deep Field (HUDF) to = 28.0 mag. The 
data are split into four sub-stacks of approximately equal exposure times. Objects exhibiting the 
signature of optical variability are selected by studying the photometric error distribution between 
the four different epochs, and selecting 622 candidates as 3. Oct outliers from the original catalog of 
4644 objects. Of these, 45 are visually confirmed as free of contamination from close neighbors or 
various types of image defects. Four lie within the positional error boxes of Chandra X-ray sources, 
and two of these are spectroscopically confirmed AGN. The photometric redshift distribution of the 
selected variable sample is compared to that of field galaxies, and we find that a constant fraction of 
~1% of all field objects show variability over the range of 0.1 < 2 < 4.5. Combined with other recent 
HUDF results, as well as those of recent state-of-the-art numerical simulations, we discuss a potential 
link between the hierarchical merging of galaxies and the growth of AGN. 

Subject headings: Galaxies: active — galaxies: formation — Active Galactic Nuclei — Supermassive 
Black Holes 



1. INTRODUCTION 

The Hubble Ultra Deep Field (HUDF: lBeckwith et all 
|2005|) is the deepest optical image of a slice of the Uni- 
verse ever observed. As such, it allows for a variety 
of different investigations into astrophysical objects and 
processes. The HUDF observations consist of 400 orbits 
with the Hubble Space Telescope (HST) over a period of 
four months in four optical bands {BVi'z') with the Ad- 
vanced Camera for Surveys (ACS), and are supplemented 
in the JH bands with an HST Legacy Program using the 
Near-Infrar ed Camera and Mu l ti-Object Sp ectrograph 
fNICMOS: IBouwens et al]l200l iThmnpson et al. 2005j). 
Given the unprecedented quality of these data, the list 
of supporting data from both space and the ground is 
constantly growing. 

Since the ACS data was observed over a period of four 
months, it provides a unique opportunity to search for 
variability in all types of objects to very faint flux levels 
(perhaps even for AB > 30 mag), such as faint stars, dis- 
tant supernovae (SNe), and weak active galactic nuclei 
(AGN). In this paper, we search for weak AGN vari- 
ability in the i'-band (F775W), because, in this filter, 
the HUDF images are deepest, and have the best tem- 
poral spacing over the four months for the variability 
search. It should be noted that at higher redshifts, the 
ACS filters sample further into the rest-frame ultraviolet. 
This is advantageous, because A GNs are known to show 
more variability in the UV fe.g.. iPaltani & Courvoisieil 
[nil). In the original Hubble D eep Field North ( HDF- 
N) and the Groth Strip Survey, iSaraiedini et all ()2003D 
and B 

araicdini (2003) performed a similar search to 
AB=27.0 mag over longer time baselines (5-7 years). 
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They searched for variability in the nuclei of the galaxies 
using small apertures, which is a different approach than 
in the present search. 

Fro m the WMAP polarization results ijKoeut et alJ 
120031) . population HI stars likely existed at z~20. These 
massive stars (> 250 Mq) are expected to produce a 
large population of massive black holes (Affc/i > 150 Mq, 
[Madau & Rees (20Q^). Since there is now good dy- 
namical evidence for the existence of supermassive 
{Mbh ^ lO^-lO^ Mq) black holes (SM BH) in the cen- 
ters of g alaxies at z ^ dK ormcn dy fc Richstond 
19951 iM agorrian. Tremainer fc Richstond Il998t 
Kormendv & Gebhardt 2001), it is important to 
understand if there is any relationship between the 
formation of the SMBHs observed at z ~ and the 
lower mass BHs at zc::; 20. A compreh e nsive review 
of SMBHs is given by iFerrarese fc FordI l|2004fl . An 
important question to address is how these SMBHs, 
as seen nearby, have grown over the course of cosmic 
time. One suggestion is that they "grow" through 
the mergers of galaxies that themselves contain less 
massive BHs, so the byproduct is a larger single galaxy 
with-eventually-a more massive BH in its center. The 
growth of the BH may then be observed via its AGN 
activity. If this scenario is true, then perhaps there 
exists an observable lin k between galaxy mergers and 
increased AGN activity fSilk fc R,eedl ll998D . Therefore, 
studying this possible link as a function of redshift could 
give insight into the growth of SMBHs and its potential 
relation to the process of galaxy assembly. 

In §2, we present the HUDF observations and sum- 
marize the essential elements of its data reduction, and 
in §3 we present the variable candidate selection. In §4 
we present the photometric redshift distribution of the 
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TABLE 1 
Observations 



Epoch 


Start 


End 


Exp. Time^ 


# of 


RJD'' 


Days Since 


No. 


Date 


Date 


(Seconds) 


Exp. 


(Days) 


Epoch 1 


1 


2003-09-24 


2003-10-10 


92340 


76 


52914.2 




2 


2003-10-10 


2003-10-29 


92340 


76 


52926.7 


12.5 


3 


2003-12-04 


2003-12-18 


89940 


76 


52985.9 


71.7 


4 


2003-12-22 


2004-01-14 


72490 


60 


53005.7 


91.5 



''There are two exposures per HST orbit 
''RJD=median Revised Juhan Date— 2.4x10^ days 



variable objects together with that of the HUDF field 
galaxies, in §5 we present the results, and in §6 we dis- 
cuss our results in terms of galaxy assembly and AGN 
growth. 

2. OBSERVATIONS 

All data use d here are from the H ubble Ultra Deep 
Field (HUDF: iBeckwith et alJ l2nol . The individual 
cosmic-ray (CR) clipp ed images and weight m aps were 
used with multidrizzle fKoekem oer et ai] l2002') to create 
eight sub-stacks of approximately equal exposure times. 
These used the same cosmic-ray maps and weight maps 
employed to create the full-depth HUDF mosaics. All 
HUDF images were drizzled onto the same output pixel 
scale (O'.'OSO per pixel) and WCS frame as the original 
HUDF. Given the time-spacing of the exposures and 
the desire to extend the study to the faintest possible 
flux- levels, the images and weight maps were combined 
in groups of two to create four epochs of observation 
for the variability study on 0.4-3.5 months timescales. 
Exposure-time weighted averages were created for all 
images, and simple addition was used to combine the 
weight maps ijBeckwith et al.ll2005l) . The exposure times 
and median observation dates are listed in Table 2] The 
four epochs chosen here have exposure times as close to 
each other as possible, so that the flux error distribu- 
tions will be as much as possible symmetric, and there- 
fore more easily modeled. As seen in Tabled this is done 
at the expense of not having the endpoints of the four 
epochs well-spaced in time. In order to be optimized 
for variability studies, future observations of this kind 
should take into account the need for both equal depth 
exposures and well separated observation dates, although 
this would further complicate the already difficult task 
of scheduling observations such as the HUDF. All magni- 
tudes are on the Ai?-system usin g the zero-points give n 
in the HUDF pubhc data release ijBeckwith et alJl2005j) . 

2.1. Catalog Generation and Photometry 

Variability was searched for by comparing the photo- 
metric catalogs from the various epochs to each other. 
Catalogs were generated using SExtractor Version 2.2.2 
(jBertin fc Arnout^ll996D with a l.Ocr detection thresh- 
old, and requiring a minimum of 15 connected pixels 
(i.e. aproximately the PSF area) at this limit above 
sky. Since we are searching for any signs of variabil- 
ity, we chose to use a liberal amount of de-blending 
{DEBLEND_MINCONT = 5x10-^). This allows for 
pieces of merging galaxies to be measured separately to 
enhance the chances of finding variable events in point- 
source components. SExtractor was run in dual image 
mode using the full-depth HUDF as the detection image. 



and utilizing the corresponding weight-maps to minimize 
the number of false detections due to edges and other 
image defects that are reflected in these weight-maps. 
This results in catalogs with 27819 measured objects, 
which still contains many over-deblended objects or edge- 
effects. Since t he HUDF was observ ed at four different 
position angles l)Beckwith et al.l2005(l . only the 15205 ob- 
jects observed in all four epochs were considered. The re- 
sult is a catalog of 12514 objects, which is 90% complete 
to i' < 30.5 mag. Since we can only measure variabil- 
ity from the individual epoch images that are one-fourth 
the full HUDF in length, the variable candidate sample 
is restricted to the 4644 objects with i' < 28.0 mag. 

The HUDF is the deepest optical image ever observed, 
and will possibly remain the deepest until the James 
Webb Space Telescope (JWST) is launched in 2011. To 
explore the limits of the HUDF depth, a few words about 
the point spread function (PSF) of the ACS images are 
needed. A comp rehensive study o f HST / ACS PSF-issues 
can be found in lHevmans et all lj2004f) . so here we only 
highlight the aspects relevant to our study. First, the 
ACS camera is not located on the optical axis of HST, 
and therefore the HUDF field is rectified by applying ge- 
ometric distortion correction polynomials. Secondly, the 
ACS/WFC PSF is known to vary with location on the 
CCD detectors, and with the time of observation due 
to "breathing" of the HST Optical Telescope Assembly 
(OTA). Since the HUDF was observed at four different 
roll angles over four months, these PSF effects can easily 
be seen by inspecting the locations of bright objects in 
an image created by dividing two images taken at dif- 
ferent roll angles. Owing to these PSF issues and the 
significantly complex ACS image registration, this "ratio- 
image" will easily show the cores of all bright objects with 
significant positive or negative flux excursions, regardless 
of whether or not they are truly variable. For this rea- 
son, we cannot use small PSF-sized ap ertures to search 
for nu cle ar variability, as was done by iSaraiedini et alJ 
()2003|) . ISaraiedini et al.l (|2003) could use the small- 
aperture method, because for the much larger WFPC2 
pixel-size and the on-axis location of the WFPC2 cam- 
era, the geometrical distortion correction and registra- 
tion effects are much smaller. Instead, we chose to use 
total magnitudes of the highly deblended objects. Even 
though our total flux apertures may encompass the whole 
galaxy, the variability necessarily must come from a small 
region (less than the 0'.'084 PSF), due to the finite hght- 
travel time across the region of variability. 

3. CANDIDATE SELECTION 

The catalogs of the four epochs were all compared to 
each other resulting in six sets of diagrams similar to the 
one shown in Fig. ^ These figures show the change in 
measured magnitudes in the SExtractor matched aper- 
tures as a function of full-depth HUDF flux. In order to 
determine which objects were true outliers (i.e., variable 
candidates), the expected error distribution for each set 
was computed as follows. For each measurement in a 
given epoch, we compute the total flux error for that i*^ 
flux measurement: 
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Fig. 1. — Magnitude difference between two HUDF epochs of 
all objects versus their total j'-band flux in AB-mag from matched 
total apertures. The ztlcr, itSo", itScr lines are shown. The black 
circles show the 222 variable candidates chosen to be \Amag\ >3o" 
outliers between any one of these these two epochs and 20 
28.0 mag. Blue points show the \Amag\ < Ha points used to 
normalize the error distribution, so that l.Ocr reflects as much as 
possible the true Gaussian l-a line. Large red points show the 
"best" 45 candidates that were chosen from all six possible epoch 
combinations, many of which were seen at > S.Ocr in 2 or more 
epoch combinations. 

where ai is the RMS per pixel in the sky background, 
A{Fi) is the number of pixels belonging to object i of a 
given flux (described in detail below) , Fi is its measured 
flux in e~ per second and G is the gain in units of seconds. 
This quantity is computed for each epoch, and for each 
of the epoch-pairs they are combined in quadrature: 



Amag = ±y^(1.0857iV)2 x ((af *)2 + (af (2) 

where N is the number of a for which the error-model is 
computed, and i,j denote the measurements at a given 
flux-level in each of the two epochs under consideration. 
Since each object at a given flux level can have a differ- 
ent area A (i.e., number of pixels), we need to assume a 
general relation between flux and area in order to opti- 
mally model the true error distribution using the above 
equations. This relation was determined iteratively for 
each pair of observations, such that 68.3% of the points 
lie within the boundaries of the upper and lower l.Ocr 
lines. We started this process by fitting the relation be- 
tween the S Extractor magnitude and ISOAREAJMAGE 
parameters as a first guess. It is then assumed that flux is 
proportional to area, and we solve iteratively for the pro- 
portionality constant to get ± l.Ocr lines that maximally 
represent a Gaussian error distribution (Fig.^. In order 
to demostrate the Gaussian nature of this error distri- 
bution at all flux levels, the Amag data are divided by 
the 1.0-cr model line, and histograms were computed for 
the resulting Amag data at various flux- levels (Fig. \^ . 
These histograms are remarkably well fit by Gaussians 
of CT ~ 1. The HUDF noise distribution is not perfectly 
Gaussian, but with 288 independent exposures in the i'- 
band, the error distribution is as close to Gaussian as 
seen in any astronomical CCD application. 

Once this Icr line is determined, we set N = 3.0 in 
Eqn. 121 and find all the objects that are at least 3. Oct 
outliers. In Fig. we show the ±1ct, ±3ct, and ±5ct 




-4 -2 2 4 -4 -2 2 4 
(Amag/Model) 

Fig. 2. — Gaussian nature of the flux error distribution at all 
flux levels. The Amag data from Fig. are divided by the best- 
fit model 1.0-cr lines. Histograms are computed from the resulting 
data for the indicated magnitude ranges. All distributions are well 
fit by Gaussians (parabolas in log space) with (t~1.0 as indicated 
in the individual panels. The almost indistinguishable dashed and 
solid lines are for the best— fit and ct = 1 Gaussians, respectively. 

lines, along with the actual data. For this particular pair 
of catalogs, there were initially 222 objects which varied 
in flux by more than 3. Oct. The choice of 3. Oct implies 
that we can expect 0.27% random contaminants. Among 
4644 objects to i' = 28.0 mag, this is about 13 interloping 
objects that are potentially classified as variable, simply 
because of the chosen 3. Oct cut-off. 

4. PHOTOMETRIC REDSHIFTS 

Photometric redshifts were computed in two ways. In 
the first method, magnitudes are computed in a fixed 
aperture with a one arcsecond radius, and objects are 
selected from the ACS z'-band. In addition to the 
HUDF ACS BVi'z' im ages, the available NICMOS JH 
l|ThomDSon et al.l2005D and VLT ISAAC K band images 
are used where available. The z'-band selection allows for 
z>5.5 galaxies to be included in the list, but the z'-band 
is not as deep as the i'-band or 1^-band images in the 
HUDF. Therefore, the primary object definition catalog 
was made in the i'-band, which introduces a bias against 
z > 5.5 objects (see Yan & Windhorst 2004). The use 
of the large apertures allows for the ground-based seeing 
K-ha,nd fiuxes to be included for more precise photo-z es- 
timates. One possible problem is that most faint galaxies 
are significantly smaller than these apertures, such that 
problems may arise for objects in crowded regions. 

In an attempt to address these issues, we also tried us- 
ing magnitudes measured within apertures defined in the 
i'-band, using the same apertures in which the variability 
was searched for [SExtractor parameter MAG-AUTO). 
The i'-band selection limits us to objects with z < 5.5. 
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This photometry is only apphed to the ACS BVi'z' 
and NICMOS JH data, which have the necessary res- 
olution to accurately measure fluxes on sub-arcsecond 
scales. The VLT K-hand data are not used here, be- 
cause they are limited by ground-based seeing (< 1" 
FWHM). The disadvantage of not using the X-band 
is the lower redshift accuracy, but the advantage is 
that the flux is measured from the same object com- 
ponent in all filters, so that crowding is less of an is- 
sue for this method. The fluxes and errors measured 
in this way are t hen input into the ph otometric red- 
shift code hvver-z |BQlzQiiella^£L_alJl2QIMl). using a suite 
of both empirical l|Coleman. Wu. fc WeedmanlllQSf]^ and 
evolutionary spectral sy nthesis (GISSEL98 update to 
iBruzual fc CharlotilT993|) templates. 

While these methods each have their own advantages 
and disadvantages, they produced the same important 
results discussed below. The difference between the 
methods are minor in the photometric redshift distri- 
bution produced, and since we discuss ratios of photo- 
metric redshift distributions in what follows, these dif- 
ferences are not relevant for the main argument below. 
The second photometric redshift determination method 
was adopted for all figures shown here. 

5. RESULTS 

5.1. Number of Variable Objects 

We find unique 622 out of 4644 variable candidates 
with i' < 28.0 mag from the six possible 2-epoch combi- 
nations. Of these, 66 are rejected as having i' > 28.0 mag 
in the final HUDF stack, leaving only 556 candidates to 
our magnitude limit. An object where just a single point 
in the light-curve is deviant would appear as an outlier 
in 3 out of 6 epoch combinations. This occurs in 25% 
of the candidates. Another 25% stand out in 2 out of 
6 pairs, and 40% stand out in 1 out of 6 pairs (usually 
indicative of a global rise or decline as a function of time 
in the light-curve). Of these 556 initial candidates, we 
find that only 45 out of 4644 show a clear sign of a com- 
pact region indicative of a point source, and are devoid 
of image defects or object splitting issues. These object 
crowding or splitting issues arise due to the extreme de- 
blending threshold (see § I2.1|l . which causes unreliable 
detections and photometric measurements of faint ob- 
jects in the wings of bright objects. Therefore, we have 
45 of 4644 objects that show the signature of AGN vari- 
ability. The 13 interlopers discussed in § O should, in 
the absence of other information, be evenly distributed 
amongst the 622 initial candidates, and therefore we ex- 
pect on the order of 1 out of 45 of our candidates to be 
a random contaminant. In order to simultaneously show 
the data from all six possible combinations of the data, 
the results are plotted in Fig. |31 showing the number of 
(7 by which each object varied in each epoch pair. The 
colored symbols are for the 45 "best" candidates. 

Another 57 objects were found that are "potentially" 
variable candidates. These are relatively isolated objects 
with reliable photometry, but show no clear sign of a 
point source. Since variability has to come from a point- 
source due to light-travel time considerations, we will 
ignore these 57 objects for now, but in §6 we will discuss 
the incompleteness resulting from the finite variability 
timcscales sampled, and from the fraction of non- variable 
or dust-obscured AGN that have likely been missed al- 
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Fig. 3. — Number of a that each object varies for all six possible 
combinations of the four different epochs. The "best" sample of 45 
variable candidates is shown with colored symbols. These objects 
were chosen to be | Amag| > 3. Oct outliers, have 20< i' <28.0 mag, 
and to be unaffected of any local image or weight map structures. 
Each object appears six times in this plot, and most candidates 
were seen at >3ct in at least two epoch pairs. Significant outliers 
not plotted in color were almost exclusively due to over-splitting 
or deblonding issues with large objects (mostly occurring in bright, 
large spiral galaxies), where the enhanced uncertainties in the local 
(object -|-sky)-subtraction introduced larger— than-Gaussian flux er- 
rors. 

together. The four-epoch light-curves for the 45 best 
candidates are shown in Fig. 01 These light-curve data 
are tabulated in Table[21 which also specifies the number 
of epoch pairs where each object was a 3. Oct outlier. For 
the 45 best candidates, 49% were discovered from a sin- 
gle pair, and 43% in two pairs. Only 5% (2 objects) were 
found in 3 pairs, which is indicative of a single deviant 
point in the 4 point light-curve. 

In Fig. 13 we show the photometric redshift distribu- 
tion for all objects with i' < 28.0 mag along with that of 
our best 45 candidates. It is clear that the distribution 
follows that of the field galaxies, i.e., there is no redshift 
where faint object variability was most prevalent. To 
make this more clear, we plot in Fig. the ratio of the 
N{z) for variables to that of field galaxies, which shows 
that this fraction is roughly constant at approximately 
1% over all redshifts probed in this study. Thi s vari- 
ability fraction is similar to the 2% found bv lSaraiedinil 
|2003) in a search for nuclear variability in the HDF-N 
and the Groth Strip survey. 

Interestingly, iStraug-hn et al.l 1)20051) show, in a com- 
panion study, that the redshift distribution for "tad- 
pole" galaxies also traces that of the HUDF field galaxies. 
They argue that these tadpole galaxies are dynamically 
unrelaxed systems, and therefore in the early stages of 
merging. The question then arises if tadpole galaxies 
and objects with point-sources that show signs of AGN 
variability are drawn from the same population. This 
question is further addressed in § |S| 

5.2. Other Tests of the Reliability of the Variable 
Sample 

5.2.1. Difference Images 

We attempted to also detect the point-source variabil- 
ity using "difference" or "ratio" images between any of 
the combinations between the two epochs. To do this, we 
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Fig. 4. — Light curves of the 45 best candidates with signs of optical AGN variabihty. The vertical axes are the change in measured 
magnitudes plotted in the sense of the average minus individual measurements. The horizontal axis is the number of days since the first 
measurement. Each panel shows the object ID number in the upper left, and the i'^g magnitude in the upper right. Plots are arranged 
in order of decreasing flux, and the combined total flux error bars are from SExtractor. The two point sources discussed in § I5.2.ijl are 
indicated by "PS." 
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ASU Photometric Redshift 

Fig. 5. — Photometric redshift distribution of all HUDF field 
galaxies (solid line) and for the "best" variable candidates (dashed 
line) multiplied by 60 X for best comparison. Photo-z's computed 
using hyperZ ^olzonella et al. 2000) and BVi'z'JH HST data 
for all galaxies with i' <28.0 mag. The redshift distribution of the 
variable objects follows that of field galaxies in general. 
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Photometric Redshift 

Fig. 6.— Percentage of HUDF objects to i' < 28.0 AB-mag 
showing variable point sources as a function of photo-z. Within the 
statistical uncertainties, about 1% of all HUDF galaxies show point 
source variability over the whole redshift range surveyed (0<z <5). 



smoothed the image by a 5x5 and 7x7 box. For each ob- 
ject, a "variance map" is computed over the four epochs 
as follows. For each pixel, we compute the nxn box- 
average for both the image and the weight map. Then 
we compute the weighted standard deviation over the 
four epochs for each pixel in the box-averaged images. 
The nxn box-size was chosen to smooth out most of 
the PSF breathing and registration issues. Even so, this 
variance image clearly had peaks in the centers of all 
bright objects. However, the variable candidates also 
tended to stand out more than their non-variable neigh- 
bors. This method worked best for the brighter objects, 
clearly verifying 13/18 candidates with < 25 mag, 
but only finding ^^50% of the full i'j^g < 28 mag sample. 
While not an optimal method of detection, this ratio- 
method provided a good way of confirming some of the 
objects whose total flux changed significantly over the 
four month period. However, the total flux method re- 
mains our primary variable candidate detection method, 



since it is most robust against ACS PSF effects. 

5.2.2. Variability Between Two Deeper Epochs 

Given the observing dates listed in Table a natu- 
ral test is to combine the first and second epochs, and 
also the third and fourth, giving only two independent 
HUDF epochs but more widely spaced in time, and with 
a somewhat higher signal-to-noise ratio for the variable 
candidate detection. However, we expect fewer variable 
candidates this way, because there are fewer epochs to 
compare, and because some of the short-term variabil- 
ity signature is necessarily smoothed over. From this 
2-epoch test, 242 candidates are chosen from the error 
distribution. Next this list is compared to the list of the 
45 best candidates from the 4-epoch search, and 12 ob- 
jects were found in both searches. If the list of the best 
and possible candidates are combined, there are 102 can- 
didates, and 30 of these are recovered in the two epoch 
test. This test gives further confidence that at least 30% 
of our top candidates are truly variable objects, although 
it does not exclude any of the other ones, since those were 
found when the shorter time-baselines were included. 

A sample of local AGN light-curves can be used in or- 
der to assess the expected completeness in temporally 
condensing our HUDF data from four to two epochs. 
The best publicly available data is that from the ex- 
tensive Inte rnational AGN Watch^ (hereafter lAGNW; 
iPeterson et al. 2002. and references therein). We used 
their five lAGNW Seyferts (all with z < 0.05), which 
have B-band light-curves which best match our observed 
HUDF i'-band data, since the median redshift of our can- 
didates is Zjried — 1-5. Since the rest-frame time sampling 
depends on the redshifts of the HUDF objects, we re- 
sampled the lAGNW light-curves for a range of redshifts 
between 0.5 < z < 5. The lAGNW data allowed for 30- 
40 independent 4-epoch light-curves to be created for 
a given redshift. These light-curves were then scaled to 
match the observed spread in Amag as a function of mag- 
nitude for our 45 HUDF candidates. A Monte Carlo test 
was then run on mock catalogs that matched the same 
magnitude distribution as our 45 HUDF candidates. On 
average half of the ^ 40 light-curves are found when our 
4-epoch variability algorithm (using same 3cr selection 
curves) is applied. The 2-epoch algorithm is then ap- 
plied to those mock candidates and the number that re- 
main variable candidates turns out to be a decreasing 
function of assumed redshift. We recover 60-70% for 
2 < 1 and 30-40% for z > A. However, this test relies 
on using 20 templates to simulate 45 objects. A visual 
inspection of the light-curves in Fig.0|shows that ideally 
a larger number of local templates is needed for a fair 
test. Nonetheless, this limited test that is possible with 
the available lAGNW data - when taken at face value - 
may be an indication that as many as one-third of our 
45 candidates are potentially false positives, although it 
is equally likely that the small number of local test tem- 
plates are not truly representative of the distant HUDF 
variability sample. Hence, we believe that the reliabil- 
ity of our faint HUDF variability sample is likely larger 
than 67%, but we will not be able to definitively say for 

^ Tabulated data from the International 
AGN Watch can be fou nd at the URL 
|http : //www . astronomy ■ ohlo- state ■ ed'u/ $ \ s im$ agnwat ch/ 
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sure without more local light-curve templates that best 
match the high-z variable objects. 

5.2.3. Other Wavelengths 

Since the UDF is in th e Chandra Deep Field-South 
fCDFS: lRosati et 311120021) . there exists deep X-ray data 
in this field, and these data are ideal for detecting AGN. 
Within the UDF, there are 16 Chandra sources (Koeke- 
moer et al. 2005, in preparation), and we detect four 
of these as variable in the optical. One of these is a 
mid-type spiral with = 21.24 mag, that is a member of 
a small group of apparently interacting galaxies. Inter- 
estingly, this object's total flux varied by less than 1%, 
but since it was one of the brightest objects in the sam- 
ple, its variability was clearly detected at > 3.0a signifi- 
cance. Two of the others are optical point sources, show- 
ing little or no visible host galaxy, with i' — 21.12 mag 
and i' — 24.79 mag. They both have measured spectro- 
scopic AGN emission-line redshifts from the GRism ACS 
Program for Extr agalactic Science survey (GRAPES; 
iPirzkal et al]l2004|) at z~3. Our BVi'z'JH photomet- 
ric redshifts put them at z = 2.6 and z — 2.9 respec- 
tively, which is rather good given that these few brighter 
AGN are the ones that are most dominated by a non- 
stellar power-law SED. The fourth candidate is barely 
detected in any of the ACS bands, and would have been 
eliminated as a spurious candidate had it not been co- 
incident with the Chandra source, and is clearly a vi- 
able object in the VLT if-band image. The detection 
of 25% of the Chandra sources as optically variable in 
the HUDF data shows that the method employed here 
is a reliable w ay of finding the AG N that are not heav- 
ily obscured. IPaoUllo et all l|200^ find that > 90% of 
the studied CDFS X-ray sources with adequate photon 
statistics show X-ray variability, and a future compari- 
son of the X-ray and optical variability amplitudes could 
provide insight into the physical mechanism that is re- 
sponsible for this variability. 

5.3. Possible Sources of Incompleteness in the Variable 
AGN Sample 

In summary, we found 45 plausibly variable objects 
with some caveats, and this number may be as high as 
102 if the "possible" variable candidates are included. 
Hence, the true variability fraction on timescales of a 
few months (rest-frame timescale few weeks to a month) 
is no more than about 1-2% of all HUDF field galaxies. 
Except for variable stars such as SNe and novae, this vari- 
ability is most likely due to an AGN given the timescales 
and distances involved. Strolger & Ricss ( 2005) found 
only one moderate rcdshift SNe in the HUDF, and thus 
SNe cannot be a significant source of contamination of 
the present sample. 

Two other possible source of incompleteness in the 
variability study must be addressed first. Non- variable 
AGN, or AGN that only vary on timescales much longer 
than 4 months, or optically obscured AGN would not 
have been detected with the current UV-optical vari- 
ability method. Sarajedini et al. (2003a) had epochs 5-7 
years apart, and found 2% of the objects to be variable. 
It is thus possible that the currently sampled times-scale 
shorter than 4 months missed about a factor 2 of all AGN 
that are only variable on longer time-scales. 



A factor of three of the brighter AGN may have been 
missed, since their UV-optical flux may be obscured by 
a dust-torus (Barthel 1989). In the AGN unification pic- 
ture, AGN cones are two-sided and their axes are ran- 
domly distributed in the sky, so that an average cone 
opening angle of co implies that a fraction l-sin(tt') of all 
AGN will point in our direction. If iu;~45° (e.g., Barthel 
1989), then every optically detected AGN (QSO) rep- 
resents 3-4 other bulge-dominated galaxies, whose AGN 
reflection cone didn't shine in our direction. Hence, their 
AGN may remain obscured by the dust-torus. Such ob- 
jects would be visible to Chandra in the X-rays or to 
Spitzer at mid-IR wavelengths, although the available 
Chandra and Spitzer data are not nearly deep enough 
to detect all HUDF objects to AB = 28 mag. (Reaching 
these depths is prohibitive in Chandra and Spitzer inte- 
gration times, and requires the next generation of X-ray 
and IR telescopes, such as Generation-X and JWST). 
At brighter flux limits, Spitzer did indeed recently find 
a significant fraction of dust-obscured AGN not seen in 
UV-optical surveys (Urrv et al. 2004). In the AGN uni- 
fication picture, the incompleteness in UV-optically se- 
lected samples due to the dust-obscuring torus would be 
as large as a factor of 3-4 ijTreister et al.ll200^ . 

6. DISCUSSION AND CONCLUSIONS 
6.1. Fraction of Variable AGN found 

Interestingly, iStraughn et al.l p005) show, in a com- 
panion study, that the redshift distribution for "tadpole" 
galaxies also traces that of the HUDF field galaxies for 
OA < z < 4.5. They argue that these tadpole galaxies 
are dynamically unrelaxed systems, and therefore in the 
early stages of merging. The question then arises if tad- 
pole galaxies and objects with point-sources that show 
signs of AGN variability are drawn from the same p opu- 
lation. At any given redshift. ISt7aughn et all f2005l ) find 
that about 6% of all HUDF galaxies appear as tadpoles, 
and they conclude that tadpole galaxies are good tracers 
of the process of galaxy assembly. 

Together with the factor of > 2 incompleteness in 
the HUDF variability sample due to the limited time- 
baseline sampled thus far, the actual fraction of weak 
AGN present in these dynamically young galaxies may 
be a factor of ^6-8 x larger than the 1% variable AGN 
fraction found through variability in the HUDF. Hence, 
perhaps as many as >6-8% of all field galaxies may host 
weak AGN, only ^ 1% of which we found here, and an- 
other > 1% could have been found if longer time-baseline 
had been available. The other factor of 3-4 x of AGN are 
likely missing because they are optically obscured, re- 
quiring the next generation of X-ray and IR telescopes. 

6.2. AGN feeding timescales compared to galaxy merger 

timescales 

We now consider if the current variability dataset can 
constrain the merging rate of SMBHs, assuming that 
SMBHs formed during hierarchical mergers of galax- 
ies containing less massive SMBHs, and if signatures of 
galaxy mergers could be related to variable AGNs. Put 
another way, can signatures of galaxy merging be related 
to our variable AGN candidates? 

A closer inspection of our data reveals that only 
1 or perhaps 2 of the variable candidates resemble 
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the tadpole galaxies of iStraughn et al.1 l)2005j) . Re- 
cent state-of-the-art hierarchical models have suggested 
that the AGN-phase only occurs in the later stages 
of a galaxy m erger, and well after it appears in the 
tadpole ph ase (Springcl. Di Matteo. fc HernauistI 120051 
Idi Matteo. Sprin gcl. & Hernquist 2005). These models 
also predict that the AGN will likely only be visible well 
o-fter (>1 Gyr) the merger induced star- formation has 
died down, implying that the fraction of dynamically 
young tadpoles that are expected to already show active 
AGN properties is relatively small. The small overlap 
between the two populations that we observe does thus 
not exclude the possibility that faint object variability 
is tracing the growth of SMBHs. However, this SMBH 
growth can only be constrained indirectly from the vari- 
ability data, and a more detailed discussion of the con- 
nection between merg ers and AGN activity is given in 
IStraughn et alJ l)2005(l and references therein. 

Given the importance of understanding the growth and 
origin of SMBHs, several lessons can be learned from this 
work in order to better design future studies of this type. 
The time-spacing of the HUDF observations -although 
as good as possible given scheduling constraints-was not 
ideal for this type of study. It is critical to re-visit 
the HUDF with HST, with the observations optimized 
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for a faint-object variability study, including covering 
time-scales of a few years. It is also essential to plan 
deeper surveys at longer wavelengths with the JWST. 
The JWST photometric and PSF stability are crucial in 
this regard, as many of our HUDF objects show signifi- 
cant variability of less than a few percent in flux. Also, a 
limiting factor in our results is the breathing of the PSF, 
along with image noise due to correcting the geometric 
distortion of ACS and slight variations of the PSF across 
the field. Therefore, JWST must have design specifi- 
cations that are capable of meeting these requirements 
in order to permit faint object variability studies to be 
done. 
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TABLE 2 
Best 45 Variability Candidates 



n5 RA (J2000) Doc (J2000) (Jab)^ Amagl" Arnag¥ AmogS^ AmagA" #^ Notcs^ 





(h m s) 


{° 




") 


(mag) 


(mag) 


(mmag) 


(mmag) 


(mmag) 


(mmag) 






10574 


03 


32 


37.92 


-27 


46 


09.1 


20.16 


20.16 


1± 


4± 


5± 


-11± 


1 




5444 


03 


32 


44.98 


-27 


47 


36.9 


20.47 


20.49 


-14± 2 


3± 2 


2± 2 


9± 2 


2 




671 


03 


32 


41.09 


-27 


48 


53.0 


20.60 


20.61 


3± 1 


12± 1 


-12± 1 


-3± 1 


1 


Two 


8870 


03 


32 


36.67 


-27 


46 


31.1 


20.94 


20.97 


-20± 3 


-3± 3 


10± 3 


13± 4 


2 


Two, CXO 


10812 


03 


32 


39.09 


-27 


46 


01.8 


21.12 


21.12 


12± 1 


-42± 1 


4± 1 


25± 1 


4 


Two, CXO, PS 


6 


03 


32 


39.54 


-27 


49 


28.4 


21.47 


21.48 


13± 1 


11± 1 


-8± 2 


-17± 2 


1 


Two 


10330 


03 


32 


37.19 


-27 


46 


08.1 


21.48 


21.51 


-11± 3 


-19± 3 


22± 3 


7± 3 


2 


Two 


9257 


03 


32 


44.28 


-27 


46 


42.3 


22.61 


22.63 


-7± 10 


14± 10 


20± 10 


-28± 12 


1 




9719 


03 


32 


36.43 


-27 


46 


32.6 


22.77 


22.82 


3± 7 


-43± 8 


32± 7 


7± 9 


1 




6489 


03 


32 


44.78 


-27 


47 


24.8 


22.81 


22.81 


25± 5 


10± 6 


-35± 6 


-1± 7 


1 




1837 


03 


32 


33.12 


-27 


48 


29.6 


23.08 


23.12 


-49± 6 


29± 6 


-7± 6 


25± 7 


2 




7474 


03 


32 


31.51 


-27 


47 


12.3 


24.21 


24.37 


97± 17 


87± 19 


-55± 20 


-147± 27 


4 


Two 


8104 


03 


32 


42.86 


-27 


47 


02.7 


24.51 


24.53 


-62± 15 


-60± 16 


95± 13 


19± 17 


2 


CXO? 


8475 


03 


32 


38.99 


-27 


46 


56.7 


24.76 


24.95 


11± 24 


-57± 28 


133± 21 


-101± 32 


2 




8145 


03 


32 


42.83 


-27 


47 


02.5 


24.79 


24.78 


-81± 6 


-88± 7 


89± 6 


67± 7 


2 


Two, CXO?, PS 


7179 


03 


32 


30.17 


-27 


47 


16.9 


24.80 


24.87 


64± 26 


50± 29 


-136± 34 


11± 33 


2 




5802 


03 


32 


37.49 


-27 


47 


31.4 


24.86 


25.46 


2± 49 


31± 51 


132± 44 


-188± 72 


1 




6603 


03 


32 


48.21 


-27 


47 


24.1 


25.02 


25.20 


75± 30 


40± 34 


-140± 37 


13± 39 


1 




8197 


03 


32 


31.83 


-27 


47 


02.9 


25.11 


25.17 


-86± 37 


19± 37 


-60± 37 


116± 38 


1 




5299 


03 


32 


36.05 


-27 


47 


37.8 


25.16 


25.18 


-56± 35 


147± 32 


1± 34 


-109± 45 


2 




6378 


03 


32 


36.20 


-27 


47 


26.2 


25.43 


25.47 


-20± 42 


119± 40 


74± 40 


-200± 62 


2 




5576 


03 


32 


42.23 


-27 


47 


33.4 


25.55 


25.56 


50± 37 


-14± 42 


87± 36 


-135± 53 


1 




9306 


03 


32 


41.09 


-27 


46 


42.4 


25.78 


25.77 


153± 33 


-252± 53 


99± 35 


-43± 49 


2 




4352 


03 


32 


33.93 


-27 


47 


49.9 


25.89 


25.94 


119± 41 


-172± 59 


-50± 49 


78± 52 


1 




6203 


03 


32 


37.00 


-27 


47 


26.3 


25.92 


26.59 


228± 73 


60± 92 


-482±141 


71±103 


3 




4094 


03 


32 


33.11 


-27 


47 


52.1 


25.96 


25.95 


-50± 34 


-247± 44 


84± 31 


169± 34 


1 


Two 


6480 


03 


32 


41.40 


-27 


47 


23.9 


25.97 


26.08 


-241± 69 


-49± 63 


104± 51 


145± 59 


1 




5692 


03 


32 


44.97 


-27 


47 


33.6 


26.03 


26.08 


191± 49 


82± 59 


-179± 71 


-137± 81 


1 


Two 


4736 


03 


32 


43.25 


-27 


47 


44.0 


26.15 


26.27 


36± 42 


81± 44 


161± 39 


-343± 73 


2 




1770 


03 


32 


35.75 


-27 


48 


31.5 


26.32 


26.30 


-24± 54 


173± 50 


-274± 72 


76± 63 


2 




8370 


03 


32 


41.97 


-27 


46 


58.1 


26.36 


26.48 


-59± 78 


186± 67 


127± 66 


-323±121 


2 




5652 


03 


32 


42.05 


-27 


47 


33.1 


26.59 


26.61 


-342± 79 


46± 60 


-105± 65 


301± 54 


2 




7882 


03 


32 


45.76 


-27 


47 


05.7 


26.73 


26.67 


115± 87 


-322±142 


238± 79 


-113±131 


2 




2511 


03 


32 


44.36 


-27 


48 


16.2 


26.76 


26.90 


-283±124 


-213±126 


66± 93 


321± 89 


2 


Two 


9286 


03 


32 


45.23 


-27 


46 


42.5 


26.90 


26.95 


-60± 94 


-22± 98 


306± 68 


-316±145 


1 




3990 


03 


32 


33.14 


-27 


47 


52.7 


26.91 


26.89 


-41± 64 


-411± 97 


127± 56 


222± 61 


1 




6251 


03 


32 


41.05 


-27 


47 


25.7 


26.98 


27.03 


-190±114 


-408±149 


217± 78 


244± 91 


2 


Two 


3670 


03 


32 


35.94 


-27 


47 


58.0 


27.04 


27.52 


199±139 


-575±309 


403±117 


-299±268 


3 




4680 


03 


32 


38.62 


-27 


47 


44.7 


27.16 


27.24 


-49±104 


20±105 


-331±139 


275± 96 


1 




7964 


03 


32 


44.27 


-27 


47 


03.6 


27.20 


27.25 


-386±156 


9±118 


346± 82 


-96±147 


1 




11336 


03 


32 


39.68 


-27 


45 


46.0 


27.22 


27.25 


-38±104 


344± 80 


-310±133 


-101±132 


1 




5936 


03 


32 


45.20 


-27 


47 


29.1 


27.49 


27.61 


-543±219 


-127±162 


213±113 


274±128 


1 




7094 


03 


32 


44.26 


-27 


47 


14.9 


27.59 


27.58 


404± 78 


104±112 


-299±153 


-404±203 


1 


Two 


2243 


03 


32 


35.93 


-27 


48 


20.5 


27.69 


27.69 


-110±168 


432±110 


133±140 


-802±396 


2 




1136 


03 


32 


37.58 


-27 


48 


46.3 


27.82 


27.90 


428±122 


142±171 


-119±211 


-793±476 


1 





Tkleasured magnitude from total UDF stack. 
''Magnitude of average flux over each of the four epochs. 

"^Difference between average magnitude and magnitude from Epoch N in miUi-mags. 

''Number of epoch pairs where variabiUty detected (max is 6). 

"^Two-also detected in two epoch test; CXO-located at position of CXO source; PS-point source (see S l5T3l . 



